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Description 

BACKGROUND OF THE INVENTION 

5 FIELD OF THE INVENTION 

The present invention relates to an ion implantation process simulation device realizing accurate interpolation and 
extrapolation of ion implantation profiles including tails in ion implantation process simulation for a semiconductor 
device, and a simulation method therefor. 

10 

DESCRIPTION QF THE RELATED ART 

A manufacturing process of semiconductor devices including LSI includes a step of implanting impurity ions as 
dopant into a semiconductor substrate by an ion implantation technique and diffusing and activating them through a 

75 thermal treatment It is well-known that an impurity distribution obtained in this step accounts for the considerable 
change in a threshold voltage Vt, an ON-state current of a transistor and other electrical characteristics parameters. In 
recent years, in particular, for reducing a semiconductor in size to achieve high density large capacity and high speed, 
designing a transistor with shallower junctions has been gaining more importance. More specifically, demanded is not 
only accurate control of a peak concentration of impurities but also accurate adjustment of position of the tail of an ion 

20 implantation profile. 

In order to meet this demand, ion implantation process simulation using a computer is adopted. This is a method 
of predicting and calculating various ion implantation profiles for use in manufacturing semiconductor devices by using 
predetermined models and algorithms on a computer. Use of ion implantation process simulation largely contributes to 
the reduction in the number of experimental steps for improving element characteristics of a semiconductor device and 

25 to the improvement in efficiency of designing work. 

In simulating ion implantation profiles, in general, ion implantation profiles are described using a Gaussian or Pear- 
son function which can obtain satisfactory enough approximation when crystals are amorphous. At that time, since ion 
implantation into a semiconductor crystal substrate reflects crystaliinity of the substrate, profiles will vary depending on 
crystal orientation of the substrate and an ion implantation angle, that is, the channeling phenomenon will occur. The 

30 channeling phenomenon tends to decrease with an increase in dose value. This is because crystals become amor- 
phous as a dose value is increased. Profiles of ion implantation into a semiconductor crystal substrate can not therefore 
be described by using a simple Gaussian or Pearson function but can be described by using a plurality of functions. 

A technique of describing a profile of ion implantation into a semiconductor crystal substrate by using a function 
and simulating the same is disclosed, for example, in "An Improved Approach to Accurately Model Shallow B and BF 2 

35 implants in Silicon" (Al F. Tasch, H. Shin, and C. Park; J. Electrochem. Soc., Vol. 136, No. 3, March 1989, pp. 810-814; 
The Electrochemical Society, Inc.). The literature recites a method of more accurately simulating experimental profiles 
of as-implanted ion impurities. The experimental profiles are results of measurement obtained by the SIMS (Secondary 
Ion Mass Spectrometry) method. SIMS is a well-known method of irradiating a surface region of a semiconductor with 
a primary ion such as an oxygen ion or a cesium ion to generate a secondary ion and conducting mass analysis of the 

40 secondary ion to measure an impurity distribution. The experimentally obtained ion implantation profile is expressed 
using two Pearson functions which represent an amorphous component and a channeling component, respectively. 
The ion implantation profile N(x) is descrtoed as a sum of the two components as shown in the following equation (1). 

N(x)=D ^f main (x)+D ^f ^(x) (1) 

45 

where a main peak amorphous component function f main (x) and a subsidiary peak channeling component function f sub 
(x) are normalized functions not dependent on a dose value, while D main represents a main peak amorphous compo- 
nent dose coefficient and D sub represents a subsidiary peak channeling component dose coefficient. A total dose D T of 
so the ion implantation prof fle N(x) is described as a sum of coefficients of the amorphous component dose and the chan- 
neling component dose as shown in the following equation (2). 



D T - jNfxJdx^D^ + D^ (2) 



The Pearson function used here is described using four kinds of moment parameters, projected range Rp, standard 
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deviation ARp, skewness y and kurtosis p. 

In the following, description will be made of a conventional ion implantation process simulation method of obtaining 
an ion implantation profile for a prescribed dose by interpolation, with reference to a flow chart of Fig. 6. With reference 
to Fig. 6, first, from table data of ion implantation prof fles for several dose values, extract moment parameters, projected 

5 range Rp, standard deviation ARp, skewness y and kurtosis 0 in two normalized functions respectively representing the 
amorphous component and the channeling component (moment parameters of a Dual Pearson function), an amor- 
phous component dose coefficient and a channeling component dose coefficient (Step 601). As a result, a Dual Pear- 
son data table is prepared. Next select parameters for doses at two points most neighboring to an arbitrary dose from 
the Dual Pearson data table (Step 602). Next, out of the selected parameters, linearly interpolate the dose-dependent 

io amorphous component dose coefficient and channeling component dose efficient with respect to dose (Step 603). 

More specifically first use SIMS experimental profile data at a plurality of dose points to extract a functional param- 
eter value and a ratio of the coefficients of an amorphous component dose at a main peak to a total dose, D main /D T . The 
functional parameter and the ratio of the coefficients of an amorphous component dose at a main peak to a total dose 
are recited for boron and BF 2 in the above-descrtoed literature. 

is Fig. 7 is a diagram showing a dose dependency as a parameter to be interpolated. In Fig. 7, the abscissa repre- 
sents a dose value, the ordinate represents a main peak amorphous component dose ratio and the polygonal line rep- 
resents a dose dependency of a main peak amorphous component dose ratio. As illustrated in Rg. 7, the number of 
dose points of the experimental data shown in the diagram are five, not so many. Ratios of the coefficients of a main 
peak amorphous component dose to the total dose at the other doses need to be interpolated for all the doses. When 

20 moment data and ratios of coefficients of a main peak amorphous component dose to the total dose (D main /D T J and 
(Dmain, i+i^T, ki) °* toe ion implantation profile at dose values D T j and D T M of two points for the interpolation or 
extrapolation with respect to an arbitrary dose value D T> a are given, linear interpolation by a conventional simulation 
method will result in describing a ratio of each coefficient of the amorphous component dose and the channel compo- 
nent dose to the total dose at the arbitrary dose value D Tt a , that is, (D,^ a/D Tf J and (D sub( a /Dj t a) , as shown in the 

25 following equation (3). 



30 



(D *a^T ( a) = (D T>1 -D TJ ) (3) 



Psub,a^T,a)= 1 -( D main^T.a) 



According to the above-descrtoed conventional ion implantation process simulation method, however, when pro- 
35 files Nj (x) and N k1 (x) at two dose points D x j and D T j+1 (D T| , < D T| i+1 ) satisfy the relationship Nj (x) < N i+n (x) , inver- 
sion might occur between the intensity of a profile N a (x) obtained by linear interpolation with respect to a dose point D Tt 
a located between the two points and the intensity of the profile N^ (x) . Shown in Fig. 8 as an example are calculation 
results obtained by a standard process simulator SUPREM-3 which is widely used at present. With reference to Rg. 8, 
the curve CP3 of the profile N a (x) obtained by linear interpolation between the profile curves CP2 and CP4 according 
40 to the above-described conventional linear interpolation method crosses with the profile curves CP1 and CP2 for larger 
doses at the tail. In other words, the profile curve CP3 is undesirably shown to have a higher concentration than those 
of the profile curves CP1 and CP2. 

Based on this result, it is possible to conduct interpolation with respect to a logarithmic value of a dose on the 
abscissa according to the following equation (4). 

45 

, n , flogD T> H^PT,a)(Pmain,^ m 

(U main,a /u T,a> = (logD T>1 -logD T ,) 1 ' 

50 (D Gub.a'D T ,a>=1 "( D maln.a'D T,a) 

While the interpolation results obtained in this case are slightly improved as compared with those of Rg. 8, undesirable 
inversion of the intensity at the tail can not be avoided. 
55 As described in the foregoing, converrtionaJ ion implantation process simulation methods have a drawback that in 
simulation of an impurity profile for junction depth, undesirable interpolation is caused at a tail of the ion implantation 
profile to prevent acquisition of accurate simulation results. 
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SUMMARY OF THE INVENTION 

An object of the present invention is to provide an ion implantation process simulation device which accurat ly inter- 
polates and extrapolates ion implantation profiles for a semiconductor device covering tails and a simulation method 
5 therefor. 

According to the first aspect of the invention, an ion implantation process simulation device which conducts ion 
implantation process simulation for a semiconductor device to conduct interpolation and extrapolation necessary for an 
ion implantation profile, comprises 

io data extracting means for extracting, from ion implantation profile data, moment parameters, projected range, 
standard deviation, skewness and kurtosis, in two normalized functions respectively representing an amorphous 
component and a channeling component an amorphous component dose coefficient and a channeling component 
dose coefficient to generate a data table, 

interpolation data obtaining means for obtaining a parameter for use in the interpolation and extrapolation of a dose 
is coefficient from said data table generated by said data extracting means. 

dose coefficient interpolating/extrapolating means for using, out of parameters obtained by said interpolation data 
obtaining means, two normalized functions respectively representing an amorphous component and a channeling 
component, and an amorphous component dose coefficient and a channeling component dose coefficient corre- 
sponding to the two functions to express a profile of ion implantation into a semiconductor crystal substrate as lin- 
(^ ; v 20 ear connection of the two functions, as well as using, at the time of descnbing said ion implantation profile for a 
dose value designated by a user, dose-independent moment parameters extracted from the ion implantation profile 
and coefficients of said linear connection dependent on dose at a plurality of dose values to interpolate and extrap- 
olate a logarithmic value of said channeling component dose coefficient with respect to logarithmic values of all 
dose values, and 

25 simulation result outputting means for outputting a simulation result obtained by the interpolation and extrapolation 
of a dose coefficient by said dose coefficient interpolating/extrapolating means. 

In the preferred construction, said dose coefficient interpolating/extrapolating means, when moment data, and an 
amorphous component dose coefficient and a channeling component dose coefficient D majn j, D 8Ub> j and ^ of 
30 the ion implantation profile at dose values D T> j and D T M of two points for conducting interpolation and extrapolation 
with respect to an arbitrary dose value Dj a are given, calculates an amorphous component dose coefficient and a 
channeling component dose coefficient D main a and D sub a at the arbitrary dose value Dy a according to the following 
equation: 

35 ^ (logP T>1 -togD Ti a )D s ubj-KlogD T,a-'og p T,i) p 

6Uba= (logD T>1 -logD Ti ) 

P main .a = P T,a " P sub, a 

40 

(^h In the preferred construction, said dose coefficient interpolating/extrapolating means, when an experimentally 

extracted moment parameter is dependent on a dose value, further interpolates and extrapolates said dose-dependent 
moment parameter with respect to logarithmic values of all doses. 

In the preferred construction, said dose coefficient interpolating/extrapolating means, when moment data, and an 
45 amorphous component dose coefficient and a channeling component dose coefficient D main t , D sub j and D sub M of 
the ion implantation profile at dose values D Tt j and D T i+1 of two points for conducting interpolation and extrapolation 
with respect to an arbitrary dose value D T a are given, calculates an amorphous component dose coefficient and a 
channeling component dose coefficient a and D sub a at the arbitrary dose value a according to the following 
equation: 

50 

D (logP T>1 -JogD j, a )D sutU+CogP T,a-'09P Tjjg sub>1 

8Ub ' a ~ <logD T>1 4ogD T>i ) 

55 D main,a ssD T,a* D sub l a 

when an experimentally extracted moment parameter is dependent on a dose value, further defines dose- 
dependent moment parameters at dose values Dj ( s and Dj ^ of two points for conducting interpolation and extrapola- 
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tion with respect to an arbitrary dose value D T a as j and ^ k1 , respectively, and calculates a moment parameter 
^ a at the arbitrary dos value Oj a according to the following equation: 

(logD Tf M -lOgD T ,a)^ m,i-K'0gD T.a-'PgP jj^j m,i+1 

5 ^m^= OogD T>1 .logD TI ) 

According to the second aspect of the invention, an ion implantation process simulation method of conducting ion 
implantation process simulation for a semiconductor device to conduct interpolation and extrapolation necessary for an 
10 ion implantation prof fle, comprising the steps of: 

extracting, from ion implantation profile data, moment parameters, projected range, standard deviation, skewness 
and kurtosis, in two normalized functions respectively representing an amorphous component and a channeling 
component, an amorphous component dose coefficient and a channeling component dose coefficient to generate 
75 a data table, 

obtaining a parameter for use in the interpolation and extrapolation of a dose coefficient from said data table gen- 
erated at said data extracting step, and 

out of parameters obtained at said interpolation data obtaining step, using two normalized functions respectively 
representing an amorphous component and a channeling component and an amorphous component dose coeffi- 

20 cient and a channeling component dose coefficient corresponding to the two functions to express a profile of ion 
implantation into a semiconductor crystal substrate as linear connection of the two functions, as well as using, at 
the time of describing said ion implantation profile for a dose value designated by a user, dose-independent 
moment parameters extracted from the ion implantation profile and coefficients of said linear connection dependent 
on dose at a plurality of dose values to interpolate and extrapolate a logarithmic value of said channeling compo- 

25 nertt dose coefficient with respect to logarithmic values of all dose values. 

In the preferred construction, at said dose coefficient interpolating/extrapolating step, when moment data, and an 
amorphous component dose coefficient and a channeling component dose coefficient D majn h D 6Ub j and D sub ^ 4 . 1 of 
the ion implantation profile at dose values D T j and Or, j+i of two points for conducting interpolation and extrapolation 
30 with respect to art arbitrary dose value Dj a are given, an amorphous component dose coefficient and a channeling 
component dose coefficient D maint a and D sub( a at the arbitrary dose value D T a are calculated according to the following 
equation: 



35 



55 



. (lOgD T >1 Wt, JP sub.rK'ogD T.a-'ogDT.i) 0 gubj+l 

8Ub ' a ° OogD Ti ^ n -logD TJ ) 



Pmain,a = D T,a ' D sub,a 

40 In the preferred construction, said dose coefficient interpolating/extrapolating step further comprises the step of, 
when an experimentally extracted moment parameter is dependent on a dose value, interpolating and extrapolating 
said dose-dependent moment parameter with respect to logarithmic values of all doses. 

In another preferred construction, at said dose coefficient interpolating/extrapolating step, when moment data, and 
an amorphous component dose coefficient and a channeling component dose coefficient D maln> j, D sub| j and D sub ^ 

45 of the ion implantation profile at dose values Dj \ and D T ^ of two points for conducting interpolation and extrapolation 
with respect to an arbitrary dose value D T a are given, an amorphous component dose coefficient and a channeling 
component dose coefficient D main a and D sub a at the arbitrary dose value D Ti a are calculated according to the following 
equation: 

50 D (logD T>1 -logD-r.a) 0 Bub,i-K lQ gP T, a - |o g D T,i)P sub,l+1 

6ub ' aD (logD T>1 -logD T(i ) 



Pmain,a = D T,a*Psub,a 

when an experimentally extracted moment parameter is dependent on a dose value, dose-dependent moment 
parameters at dose values D T j and D T> k1 of two points for conducting interpolation and extrapolation with respect to 
an arbitrary dose value D T a are defined as ^ \ and Mm ki . respectively, and a moment parameter ^ a at the arbitrary 
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dose value D T a is calculated according to the following equation: 

(logD T M -logD T a )ji ml +(logD Ta -logD T ,) u 
= (logD T>1 -logD Tl ) 

5 

According to the third aspect of the invention, an ion implantation process simulation method of conducting ion 
implantation process simulation for a semiconductor device to conduct interpolation and extrapolation necessary for an 
ion implantation profile, comprising the steps of: 

10 

using two normalized functions respectively representing an amorphous component and a channeling component, 
and an amorphous component dose coefficient and a channeling component dose coefficient corresponding to the 
two functions to express a profile of ion implantation into a semiconductor crystal substrate as linear connection of 
the two functions, and 

is using dose-independent moment parameters extracted from the ion implantation profile and coefficients of said lin- 
ear connection dependent on dose at a plurality of dose values to interpolate and extrapolate a logarithmic value 
of said channeling component dose coefficient with respect to logarithmic values of all dose values. 

In the preferred construction, at said step of interpolating and extrapolating a logarithmic value of said channeling 
20 component dose coefficient with respect to logarithmic values of all dose values, with respect to an arbitrary dose value 
as a target of interpolation and extrapolation, a channeling component dose coefficient for the arbitrary dose value is 
interpolated and extrapolated by using said amorphous component dose coefficients and channeling component dose 
coefficients at dose values of two points most neighboring to the dose value. 

According to another aspect of the invention, an ion implantation process simulation method of conducting ion 
25 implantation process simulation for a semiconductor device to conduct interpolation and extrapolation necessary for an 
ion implantation profile, comprising the steps of: 

using two normalized functions respectively representing an amorphous component and a channeling component, 
and an amorphous component dose coefficient and a channeling component dose coefficient corresponding to the 
30 two functions to express a profile of ion implantation into a semiconductor crystal substrate as linear connection of 
said two functions, and 

using moment parameters extracted from the ion implantation profile and linear connection coefficients at a plurality 
of dose values to, when an experimentally extracted moment parameter is dependent on a dose value, interpolate 
and extrapolate said dose-dependent moment parameter with respect to logarithmic values of all dose values. 

35 

In the preferred construction, at said step of interpolating and extrapolating a logarithmic value of said channeling 
component dose coefficient with respect to logarithmic values of all dose values, with respect to an arbitrary dose value 
as a target of interpolation and extrapolation, dose-dependent said moment parameters at dose values of two points 
most neighboring to the dose value are used to interpolate and extrapolate a moment parameter for the dose value. 
40 According to a still further aspect of the invention, a computer readable memory having a control program for con- 
trolling an ion implantation process simulation device to conduct ion implantation process simulation for a semiconduc- 
tor device and conduct interpolation and extrapolation necessary for an ion implantation profile, said control program 
comprising the steps of: 

45 from ion implantation profile data, extracting moment parameters, projected range, standard deviation, skewness 
and kurtosis, in two normalized functions respectively representing an amorphous component and a channeling 
component, an amorphous component dose coefficient and a channeling component dose coefficient to generate 
a data table, 

obtaining a parameter for use in the interpolation and extrapolation of a dose coefficient from said data table gen- 

so erated at said data extracting step, and 

out of parameters obtained at said interpolation data obtaining step, using two normalized functions respectively 
representing an amorphous component and a channeling component and an amorphous component dose coeffi- 
cient and a channeling component dose coefficient corresponding to the two functions to express a profile of ion 
implantation into a semiconductor crystal substrate as linear connection of the two functions, as well as using, at 

55 the time of describing said ion implantation profile for a dose value designated by a user, dose-independent 
moment parameters extracted from the ion implantation profile and coefficients of said linear connection dependent 
on dose at a plurality of dose values to interpolate and extrapolate a logarithmic value of said channeling compo- 
nent dose coefficient with respect to logarithmic values of all dose values. 
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Other objects, features and advantages of the present invention will become clear from the detailed description 
given herebelow. 

pRIFF DESCRIPTION OF THE DRAWINGS 

5 

The present invention will be understood more fully from the detailed description given herebelow and from the 
accompanying drawings of the preferred embodiment of the invention, which, however, should not be taken to be limi- 
tative to the invention, but are for explanation and understanding only. 

In the drawings: 

10 

Fig. 1 is a block diagram showing structure of an ion implantation simulation device according to one embodiment 
of the present invention. 

Fig. 2 is a flow chart showing operation of the ion implantation simulation device according to the present embodi- 
ment 

is Fig. 3 is a diagram of profiles showing the state of interpolation of ion implantation profiles according to the present 
embodiment 

Fig. 4 is a diagram showing a channeling component dose coefficient according to the present embodiment as 
compared with a conventional interpolation result 

Fig. 5 is a flow chart showing another operation of the ion implantation simulation device according to the present 
20 embodiment 

Fig. 6 is a flow chart showing a procedure of conventional ion implantation simulation. 

Rg. 7 is a diagram showing dose dependencies of parameters to be interpolated in a conventional ion implantation 
process simulation method. 

Rg. 8 is a diagram of profiles showing the state of interpolation of ion implantation profiles according to the conven- 
25 tional ion implantation process simulation method. 

DESCRIPTION OF THE PREFERRED EMBODIMENT 

The prefened embodiment of the present invention will be discussed hereinafter in detail with reference to the 

30 accompanying drawings. In the following description, numerous specific details are set forth in order to provide a thor- 
ough understanding of the present invention, ft will be obvious, however, to those skilled in the art that the present 
invention may be practiced without these specific details. In other instance, well-known structures are not shown in 
detail in order to unnecessary obscure the present invention. 

Fig. 1 is a block diagram showing structure of an ion implantation process simulation device according to one 

35 embodiment of the present invention. Trie ion implantation process simulation device of the present embodiment is 
implemented by a computer system such as a work station or a personal computer, and functions of the present embod- 
iment which will be described in the following are executed by a control program which controls the computer system. 
The control program is provided as storage in a storage medium such as a magnetic disk or a semiconductor memory 
and load of the control program into a processing device of the computer system realizes the ion implantation process 

40 simulation device of the present embodiment. 

With reference to Rg. 1 , the ion implantation process simulation device of the present invention includes a Dual 
Pearson data extracting unit 10 for extracting desired Dual Pearson data from ion implantation profile data 50 to gener- 
ate a Dual Pearson data table 1 1. a Dual Pearson data for interpolation obtaining unit 20 for obtaining Dual Pearson 
parameters for use in interpolating and extrapolating dose coefficients from the Dual Pearson data table 1 1 generated 

45 by the Dual Pearson data extracting unit 10, a dose coefficient interpolating/extrapolating unit 30 for interpolating or 
extrapolating dose coefficients by using a parameter obtained by the Dual Pearson data for interpolation obtaining unit 
20, and a simulation result outputting unit 40 for outputting a simulation result obtained by interpolation of dose coeffi- 
cients. In Rg. 1 , illustration is made only of a characteristic part of the structure of the present embodiment and that of 
the remaining common part is omitted. 

so Prior to the execution of simulation, the Dual Pearson data extracting unit 10 receives input of the ion implantation 
profile data 50 for various dose values to generate the Dual Pearson data table 1 1 in advance. Stored in the Dual Pear- 
son data table 1 1 are moment parameters, projected range Rp, standard deviation ARp, skewness y and kurtosis 0 in 
two normalized functions respectively representing an amorphous component and a channeling component (moment 
parameters of a Dual Pearson function), an amorphous component dose coefficient and a channeling component dose 

55 coefficient 

The Dual Pearson data for interpolation obtaining unit 20 obtains parameters at dose values of two points most 
neighboring to an arbitrary dose designated by a designation means not shown (hereinafter referred to as a Dual Pear- 
son parameter) from the Dual Pearson data table 1 1 . 
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The dose coefficient int rpolating/extrapolating unit 30, when out of the Dual Pearson parameters obtained by the 
Dual Pearson data for interpolation obtaining unit 20, an amorphous component dose coefficient and a channeling com- 
ponent dose coefficient are dependent on a dose value, interpolates or extrapolates the amorphous component dose 
coefficient and the channeling component dose coefficient wfth respect to all the dose values in the following manner. 

5 More specifically, the dose coefficient interpolating/extrapolating unit 30 first uses two normalized functions respec- 

tively representing the amorphous component and the channeling component, and an amorphous component dose 
coefficient and a channeling component dose coefficient corresponding to the two functions to express a profile of ion 
implantation into a semiconductor crystal substrate as linear connection of the two functions. Then, at the time of 
describing an ion implantation profile for a dose value designated by a user, the unit 30 uses dose-independent moment 

10 parameters extracted from the ion implantation profile and coefficients of the linear connection dependent on dose at a 
plurality of dose values to interpolate or extrapolate a logarithmic value of a channeling component dose coefficient with 
respect to logarithmic values of all the dose values. 

The foregoing operation will be expressed by an equation. When moment data, and an amorphous component 
dose coefficient and a channeling component dose coefficient D^n it is D subj j and D sub i+1 of the ion implantation pro- 

15 file at dose values Dj j and Dj ^ of two points for the interpolation or extrapolation with respect to an arbitrary dose 
value Dx t a are given, an amorphous component dose coefficient and a channeling component dose coefficient D maini 
a and D S u b a at the arbitrary dose value Dj, a is calculated according to the following equation (5). 

n (logD w -logD T, a )P subJ+QogP T.a-'PgP T,i) D sub.i+1 

20 u sub.a= (logD T>1 -logD Ti ) w 

D main .a = D T.a " D sub,a 

25 The simulation result outputting unit 40 generates such a profile diagram as shown in Fig. 3 by using a simulation 
result obtained by the interpolation and extrapolation by the dose coefficient interpolating/extrapolating unit 30 and dis- 
plays the same on a display device or prints out the same. 

Description will be next made of operation of the present embodiment with reference to the flow chart shown in Fig. 
2. The Dual Pearson data extracting unft 10 first extracts moment parameters, an amorphous component dose coeffh 

30 cient and a channeling component dose coefficient of a Dual Pearson function from the ion implantation profile data 50 
for several dose values to generate the Dual Pearson data table 1 1 (Step 201). 

Next the Dual Pearson data for interpolation obtaining unit 20 obtains Dual Pearson parameters for dose values at 
two points most neighboring to an arbitrary dose designated by a designation means not shown from the Dual Pearson 
data table 11 (Step 202). 

35 Next, the dose coefficient interpolating/extrapolating unit 30 determines whether out of the Dual Pearson parame- 
ters obtained by the Dual Pearson data for interpolation obtaining unit 20, an amorphous component dose coefficient 
and a channeling component dose coefficient are dependent on a dose value or not and when the determination is 
made that they are dose dependent interpolates or extrapolates the amorphous component dose coefficient and the 
channeling component dose coefficient with respect to the dose value in question in the above-described manner (Step 

40 203). 

Based on the calculation results of ion implantation profiles thus obtained by interpolation and extrapolation of dose 
coefficients, the simulation result outputting unit 40 generates and outputs a profile diagram. Fig. 3 shows calculation 
results of ion implantation profiles for BF 2 according to the present embodiment 

With reference to Fig. 3, accurate profiles are obtained in which a curve CP3 of a profile N a (x) obtained by the inter- 
45 polation with respect to a dose D T> a between two kinds of doses of profile curves CP2 and CP4 crosses with neither of 
the profile curves CP1 and CP2 with larger doses at any part of the entire curve including the tail. 

Fig. 4 is a diagram showing a channeling component dose coefficient as compared with that of conventional simu- 
lation. In Fig. 4, the solid line denotes a relationship between a logarithmic value of a channeling component dose coef- 
ficient and a logarithmic value of all dose values in the simulation according to the present embodiment, while the 
so broken line denotes a relationship between a logarithmic value of a channeling component dose coefficient and a log- 
arithmic value of all dose values in conventional simulation. With reference to Fig. 4, a rugged part appears in the con- 
ventional simulation as indicated by the broken line, while in the simulation according to the present embodiment the 
relationship shows monotonic increase as indicated by the solid line. 

As described in the foregoing, the present embodiment enables accurate interpolation of ion implantation profiles 
55 including tails which is demanded in simulation for a junction depth crucial to the manufacture of semiconductor 
devices. 

Description will be next made of an ion implantation process simulation device according to another embodiment 
of the present invention. Structure of the present embodiment is the same as that of the first embodiment shown in Fig. 
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1 and therefore no illustration is made thereof. 

In the present embodiment the dose coefficient interpolating/extrapolating unit 30. similarly to the first embodi- 
ment, interpolates or extrapolates an amorphous component dose coeffici nt and a channeling component dose coef- 
ficient with respect to a dose value, and also determines whether out of Dual Pearson parameters, moment parameters 

5 of a Dual Pearson function for dose values at two points most neighboring to a dose value to be calculated are dose 
dependent and when the determination is made that they are dose dependent interpolates or extrapolates the moment 
parameters of the Dual Pearson function with respect to a dose in a manner set forth in the following. 

More specifically, when the extracted moment parameters are dose dependent the unit 30 interpolates or extrap- 
olates the above-described dose-dependent moment parameters with respect to a logarithmic value of all doses. 

io This operation will be expressed by an equation. Dose-dependent moment parameters of a Dual Pearson function 
at dose values Dj j and D T i+1 of two points for the interpolation or extrapolation with respect to an arbitrary dose value 
Dj t a are defined as ^i m> \ and ^ M , respectively. A moment parameter ^ a of the Dual Pearson function at the arbi- 
trary dose value Dj a will be accordingly expressed by the following equation (6). 

15 (lOgD t, kl - |0 9° T,a) H m,i-K'OgP T.a-^P T, j) * m,>+1 ,~ 

*V* = (logD T>1 -togD Tl ) w 



a 



Fig. 5 is a flow chart showing another example of operation of the present embodiment With reference to Fig. 5, 

20 operation at Steps 501 to 503 by the dose coefficient interpolating/extrapolating unit 30 is the same as that at Steps 201 
to 203 of the operation example shown in Fig. 3. Thereafter, the dose coefficient interpolating/extrapolating unit 30 
determines whether out of the Dual Pearson parameters, moment parameters of a Dual Pearson function for dose val- 
ues at two points most neighboring to a dose value to be calculated are dose dependent and when they are dependent, 
interpolates or extrapolates the moment parameters of the Dual Pearson function with respect to a dose in the following 

25 manner (Steps 504 and 505). 

When arsenic ions are experimentally implanted into a semiconductor crystal substrate, for example, obtained as 
skewness y is a dose-dependent parameter. In this case, according to a conventional dose interpolation method using 
linear interpolation, profile curves might cross with each other at tails in some cases, while the present invention ena- 
bles accurate profiles to be obtained in which none of prof fle curves crosses with each other at tails. 

30 As described in the foregoing, by interpolating and extrapolating a logarithmic value of a channeling component 
dose coefficient and a dose-dependent moment parameter with respect to a logarithmic value of all doses, the present 
invention enables accurate interpolation and extrapolation of ion implantation profiles including tails. 

Although the invention has been illustrated and described with respect to exemplary embodiment thereof, it should 
be understood by those skilled in the art that the foregoing and various other changes, omissions and additions may be 

35 made therein and thereto, without departing from the spirit and scope of the present invention. Therefore, the present 
invention should not be understood as limited to the specific embodiment set out above but to include all possible 
embodiments which can be embodies within a scope encompassed and equivalents thereof with respect to the feature 
set out in the appended claims. 

40 Claims 

1 . An ion implantation process simulation device which conducts ion implantation process simulation for a semicon- 
ductor device to conduct interpolation and extrapolation necessary for an ion implantation profile, comprising: 

45 data extracting means (10) for extracting, from ion implantation profile data, moment parameters, projected 

range, standard deviation, skewness and kurtosis, in two normalized functions respectively representing an 
amorphous component and a channeling component an amorphous component dose coefficient and a chan- 
neling component dose coefficient to generate a data table (1 1); 

interpolation data obtaining means (20) for obtaining a parameter for use in the interpolation and extrapolation 
so of a dose coefficient from said data table (1 1 ) generated by said data extracting means (10); 

dose coefficient interpolating/extrapolating means (30) for using, out of parameters obtained by said interpola- 
tion data obtaining means (20), two normalized functions respectively representing an amorphous component 
and a channeling component, and an amorphous component dose coefficient and a channeling component 
dose coefficient corresponding to the two functions to express a profile of ion implantation into a semiconductor 
55 crystal substrate as linear connection of the two functions, as well as using, at the time of describing said ion 

implantation profile for a dose value designated by a user, dose-independent moment parameters extracted 
from the ion implantation profile and coefficients of said linear connection dependent on dose at a plurality of 
dose values to interpolate and extrapolate a logarithmic value of said channeling component dose coefficient 
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with respect to logarithmic values of all dose values; and 

simulation result outputting means (40) for outputting a simulation result obtained by the interpolation and 
extrapolation of a dose coefficient by said dose coefficient interpolating/extrapolating means (30). 

5 2. The ion implantation process simulation device as set forth in claim 1 , wherein 

said dose coefficient interpolating/extrapolating means (30), when moment data, and an amorphous compo- 
nent dose coefficient and a channeling component dose coefficient D main jp D sub j and D sub of the ion 
implantation profile at dose values Dj - t and Dj • H . 1 of two points for conducting interpolation and extrapolation 
10 with respect to an arbitrary dose value Dj a are given, calculates an amorphous component dose coefficient 

and a channeling component dose coefficient a and a at the arbitrary dose value Dj, a according to 
the following equation: 

n (logP T>1 -logD Ta )D Mb.rKtopP T.a-'ogP t,i)P m*M 
75 8Ub ' a = (logD T h . 1 -logD T , 

D maln^ = D T,a " D sub, a 

20 3. The ion implantation process simulation device as set forth in claim 1 , wherein 

said dose coefficient interpolating/extrapolating means (30), when an experimentally extracted moment 
parameter is dependent on a dose value, further interpolates and extrapolates said dose-dependent moment 
parameter with respect to logarithmic values of all doses. 

4. The ion implantation process simulation device as set forth in claim 1 , wherein 



25 



said dose coefficient interpolating/extrapolating means (30), when moment data, and an amorphous compo- 
nent dose coefficient and a channeling component dose coefficient 

Dmain, \> Psub i Dsub, i+1 Of the ion 

30 implantation profile at dose values Dr j and Dj M of two points for conducting interpolation and extrapolation 

with respect to an arbitrary dose value Dj a are given, calculates an amorphous component dose coefficient 
and a channeling component dose coefficient D^^ a and a at the arbitrary dose value D T a according to 
the following equation: 

35 (logD T>1 -fogD T a )D ^,-KlogD Ta -logD T j)D 

sub, i+1 

6Ub - a= OogD TW -logD Xl ) 



40 



45 



P main,a ~ P T.a " P sub, a 

when an experimentally extracted moment parameter is dependent on a dose value, further defines 
dose-dependent moment parameters at dose values Dj j and D T M of two points for conducting interpolation 
and extrapolation with respect to an arbitrary dose value Oj t a as u Tn> j and M , respectively, and calculates 
a moment parameter a at the arbitrary dose value Dj a according to the following equation: 



(logD -^OgP T,a)H m.i'K'ogDT.a^OgDT.i)^ m,i+1 

^ m<a= (logD Tkr logD Ti ) 

so 5. An ion implantation process simulation method of conducting ion implantation process simulation for a semicon- 
ductor device to conduct interpolation and extrapolation necessary for an ion implantation profile, comprising the 
steps of: 

extracting, from ion implantation profile data, moment parameters, projected range, standard deviation, skew- 
55 ness and kurtosis, in two normalized functions respectively representing an amorphous component and a 

channeling component, an amorphous component dose coefficient and a channeling component dose coeffi- 
cient to generate a data table; 

obtaining a parameter for use in the interpolation and extrapolation of a dose coefficient from said data table 
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generated at said data extracting step; and 

out of parameters obtained at said interpolation data obtaining step, using two normalized functions respec- 
tively representing an amorphous component and a channeling component and an amorphous component 
dose coefficient and a channeling component dose coefficient corresponding to the two functions to express a 
profile of ion implantation into a semiconductor crystal substrate as linear connection of the two functions, as 
well as using, at the time of describing said ion implantation profile for a dose value designated by a user, dose- 
independent moment parameters extracted from the ion implantation profile and coefficients of said linear con- 
nection dependent on dose at a plurality of dose values to interpolate and extrapolate a logarithmic value of 
said channeling component dose coefficient with respect to logarithmic values of all dose values. 

The ion implantation process simulation method as set forth in claim 5, wherein 

at said dose coefficient interpolating/extrapolating step, when moment data, and an amorphous component 
dose coefficient and a channeling component dose coefficient v D sub j and D 8Ub }+1 of the ion implanta- 
tion profile at dose values D T j and D T k1 of two points for conducting interpolation and extrapolation with 
respect to an arbitrary dose value D T( a are given, an amorphous component dose coefficient and a channeling 
component dose coefficient D mairii a and D sub( a at the arbitrary dose value Oj t a are calculated according to the 
following equation: 

n (togD T m -logD T a )D g^i+QogD T , a -\ogD T j)P subj+i 
8UM " (logD T M -logD T s ) 



D main,a = D T.a ' D sub,a 

The ion implantation process simulation method as set forth in claim 5, wherein 

said dose coefficient interpoiating/extrapolating step further comprises the step of, when an experimentally 
extracted moment parameter is dependent on a dose value, interpolating and extrapolating said dose-depend- 
ent moment parameter with respect to logarithmic values of all doses. 

The ion implantation process simulation method as set forth in claim 5, wherein 

at said dose coefficient imerpolating/extrapolating step, when moment data, and an amorphous component 
dose coefficient and a channeling component dose coefficient D majni \, D subt j and D 6Ubf j + i of the ion implanta- 
tion profile at dose values D T( j and Dr ( M o* two points for conducting interpolation and extrapolation with 
respect to an arbitrary dose value D T a are given, an amorphous component dose coefficient and a channeling 
component dose coefficient D maini a and D sub a at the arbitrary dose value Dj f a are calculated according to the 
following equation: 

(logD T M -logD T a )D su^-KlogD T ,a- to 9 D tj) d sub.ki 
8Uba= (logD T>r logD TI ) 



D main,a IsD T.a- D sub,a 

when an experimentally extracted moment parameter is dependent on a dose value, dose-dependent 
moment parameters at dose values Dj t ; and D T| j+1 of two points for conducting interpolation and extrapolation 
with respect to an arbitrary dose value Dj, a are defined as j and |x m> ^, respectively, and a moment param- 
eter a at the arbitrary dose value Dj a is calculated according to the following equation: 

('POP T, k1 -'ogD T.a)^ mJ+flogD T,a- |o g D T,i)^m,M-1 

^.a = (logD T>1 -logD Tj ) 



An ton implantation process simulation method of conducting ion implantation process simulation for a semicon- 
ductor device to conduct interpolation and extrapolation necessary for an ion implantation profile, comprising the 
steps of: 
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using two normalized functions respectively representing an amorphous component and a channeling compo- 
nent and an amorphous component dose coefficient and a channeling component dose coefficient corre- 
sponding to the two functions to express a profile of ion implantation into a semiconductor crystal substrate as 
linear connection of the two functions; and 

using dose-independent moment parameters extracted from the ion implantation profile and coefficients of 
said linear connection dependent on dose at a plurality of dose values to interpolate and extrapolate a logarith- 
mic value of said channeling component dose coefficient with respect to logarithmic values of all dose values. 

1 0. The implantation process simulation method as set forth in claim 9, wherein 

at said step of interpolating and extrapolating a logarithmic value of said channeling component dose coeffi- 
cient with respect to logarithmic values of all dose values, with respect to an arbitrary dose value as a target of 
interpolation and extrapolation, a channeling component dose coefficient for the arbitrary dose value is inter- 
polated and extrapolated by using said amorphous component dose coefficients and channeling component 
dose coefficients at dose values of two points most neighboring to the dose value. 

11. An ion implantation process simulation method of conducting ion implantation process simulation for a semicon- 
ductor device to conduct interpolation and extrapolation necessary for an ion implantation profile, comprising the 
steps of: 

using two normalized functions respectively representing an amorphous component and a channeling compo- 
nent and an amorphous component dose coefficient and a channeling component dose coefficient corre- 
sponding to the two functions to express a profile of ion implantation into a semiconductor crystal substrate as 
linear connection of said two functions; and 

using moment parameters extracted from the ion implantation profile and linear connection coefficients at a 
plurality of dose values to, when an experimentally extracted moment parameter is dependent on a dose value, 
interpolate and extrapolate said dose-dependent moment parameter with respect to logarithmic values of all 
dose values. 

1 2. The implantation process simulation method as set forth in daim 1 1 , wherein 

at said step of interpolating and extrapolating a logarithmic value of said channeling component dose coeffi- 
cient with respect to logarithmic values of all dose values, with respect to an arbitrary dose value as a target of 
interpolation and extrapolation, dose-dependent said moment parameters at dose values of two points most 
neighboring to the dose value are used to interpolate and extrapolate a moment parameter for the dose value. 
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